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Abstract 



A new approach to the implementation of a quantum computer by high-resolution nuclear 
magnètic resonance (NMR) is described. The key feature is that two or more line-selective 
radio-frequency pulses are applied simultaneously. A three-qubit quantum computer has 
been investigated using the 400 MHz NMR spectrum of the three coupled protons in 2,S- 
dibromopropanoic acid. It has been employed to implement the Deutsch-Jozsa algorithm 
for distinguishing between constant and balanced functions. The extension to systems 
containing more coupled spins is straightforward and does not require a more protracted 
experiment. 
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1 Introduction 



While there has long been theoretical interest in the notion of a quantum computer, it 
was the series of recent results leading to the remarkable algorithm of Shor |j^] for finding 
prime factors in polynomial time which led to the recent explosion of interest in the sub- 
ject. These theoretical results have led many groups to try to realise a quantum computer 
experiment ally. Nuclear magnètic resonance offers a particularly attractive implementa- 
tion of quantum computers because nuclear spins are relatively weakly coupled to the 
environment, and there is a long history of development of experimental techniques for 
manipulating the spins using radio frequency pulses. 

A number of groups have already demonstrated the use of NMR computers |§, ^, 
|, I, g, |, |, |Ï0|. One of the key challen ges is to try to increase the size of the system 
used. Previous work on implementing quantum algorithms has focussed on two algorithms 



in particular, the Deutsch-Jozsa |11] algorithm for distinguishing between balanced and 
constant functions and Grover's algorithm [ |12| | for searching a database. Previous work 
on both of these algorithms has used NMR computers with two qubits. In this paper 
we take the study further by implementing the Deutsch-Jozsa algorithm for a system of 
three qubits. A particularly notable feature of the experiments we describe is the use of 
simultaneous line-selective pulses to implement the key stage of the algorithm, quantum 
gates which are closely related to the controlled-controlled-not gate. 

The Deutsch-Jozsa algorithm which we will implement is to distinguish between two classes 
of two-bit binary functions: 

/:{0,1}X{0,1}^{0,1}. (1) 

The two classes are the constant functions, in which all input vàlues get mapped to the 
same output value, and the balanced functions in which exactly two of the inputs get 
mapped to 0. The eight balanced or constant functions are given in Table |^. 

The point of the Deutsch-Jozsa algorithm is that it is possible to decide whether a function 
is constant or balanced with only one evaluation of the function /. 

The theoretical steps of the quantum algorithm are as follows: 

[1] Preparation: Prepare the system in the (pure) state ipi = |0)|0)|0). 

[2] Excitation: Perform rotations of the spins about the y-axis so that the state becomes 
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V'2 = (lo) + |i))(lo) + |i))(lo)-|i)). 

[3] Evaluation: This is done by implementing the unitary transformation 

\i)\j)\k)^\i)\j)\k + f{i,j)), 
where the addition is performed modulo two. The three qubits are now in the state 
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-i/(^·^·)|i)b·)(|o)-|i)). 



(3) 



For example in the case of the function f^, the state is V's = — (|0) — 11))(|0) + 11))(|0) — |1)). 
The function /4 is implemented by applying the unitary operator 
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to the state. We note that this may be written as 



( 2a^ \ 

20":,; 
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(5) 



where Ox is the Pauh matrix normahzed so that tr(í7^) = ^ and E is the 2x2 identity 
matrix. For convenience we will denote such block diagonal matrices by the symbol A, so 
that we write 



t/4 = A(2c7^.,2cj^,S,S). 



(6) 



The complete list of unitary operators corresponding to the eight balanced or constant 
functions is given in Table |. 

[4] Observation: By rotating back by the inverse of the transformation applied in the ex- 
citation stage, it may be noted that the outputs from the unitary operators corresponding 
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to the constant functions /i and /2 have a component proportional to |0)|0)|0), whereas 
the outputs from the balanced functions fs-.-fs are orthogonal to this vector so that the 
constant and balanced functions may be distinguished from each other with probabihty 
one by a von Neumann measurement. 

The ensemble nature of an NMR quantum computer means that the implementation of 
the algorithm differs somewhat from the theoretical version. In particular the prepara- 
tion stage differs since the states of the system are not pure and the observation stage 
does not use a von Neumann measurement but measures the amplitudes of spectral lines. 
Nonetheless the key goal of the algorithm remains to determine whether the unitary op- 
erator which acts on the system in the third step corresponds to a constant or balanced 
function. 

One way to proceed would be to follow and produce a pseudo-pure state in the prepa- 
ration stage. In terms of product operators, the state corresponding to the pure state 
V'l = |0)|0)|0) is written 

p{i;i) = I, + S, + + 2I,S, + 2I,R, + 2S,R^ + 4I^S,R^. (7) 

I refers to the first spin, S, the second and R, the third. 

This state is excited to p{'ip2) = Ix + — Rx + "^^x^x ~ 2IxRx ~ "^^x^x ~ ^I^S^Rx in 
stage [2]. 

In the evaluation stage, the state to which the spins evolve depends on which function 
is being implemented. For example the unitary operator corresponding to produces 
p{ip3) = —Ix + Sx — 'Rx — '^ïxSx + '^ïx^x — '^Sx^x + ^ïxSx'í^x- The full list of output states 
is given in Table |3|. 

It should be noted that, of the observable terms (i.e. those terms linear in I^^, and 
Rj^), the term in R^; always has the same phase, but the balanced functions have altered 
signs of I^. or S^, or both. Thus if one observes that the Ix or S^; (or both) quartets are 
inverted one knows that the function is balanced. 

We note however that the same goal can be achieved by starting with thermal rather than 
pure initial states. This is because, as we will show below, similar effects are observed from 
the outputs starting with thermal initial states as were visible starting from pure initial 
states. This is not the first time that it has been noted that in NMR quantum computers. 
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thermal initial states are sufíicient to implement the algorithms of interest . 

Thus in the NMR implementation that we will use, the theoretical steps [1] to [4] are 
replaced with 

[1*] Preparation: One starts with the thermal initial state 

I, + S, + (8) 

[2*] Excitation: Apply a hard 7r/2 pulse along the y-axis to arrive at 

+ + (9) 

[3*] Evaluation: Now evolve the system with one of the unitary operators given in Table 
|2[ This is achieved by using simultaneous line selective pulses (see below). 

For example under the state evolves to 

2I^R^ + + R^. (10) 

The list of states to which each of I^;, S^; and R^^ evolve is given in Table ^. 

However (see below) the line selective pulses produce evolution by a unitary operator 
which is close to that required but differs by a controlled phase shift. For example, in the 
case of /4, the line selective pulse produces the unitary transformation 

A{2ia^,2ia^,E,E), (11) 

whereas the unitary given in Table ^ is 

U4 = Ai2a,,2a,,E,E). (12) 

The relation between these two matrices is 

Ui = A{2aa:,2a^,E, E) = A{2ia^, 2ia^,E, E) x A{-iE, -iE, E, E). (13) 

The second matrix on the right hand side of this equation is a z rotation on the first spin 
by the angle 7r/2. Thus if one wants to implement [/4, it would be necessary to follow the 
line-selective pulse by a phase shift. One finds that similar phase shifts are required for 
all functions except /i and f2- 
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[4*] Observation: Under evolution by the unitary operators corresponding to any of the 
balanced functions, either the / response or the S response (or both) disappears. Had we 
started with a pure initial state the equivalent hne would have been inverted. 

We note that the disappearance or otherwise of the I or S response is not affected by the 
final phase shift. This is because the state 

+ + (14) 

still evolves to states in which the same line disappears even if this last phase shift is not 
implemented. This may be appreciated by looking at the product operators to which the 
state evolves, as given in Table ^. 



2 Experimental Realization 

One possible way to implement the evaluation stage of the algorithm would be to make 



use of the fact [13| that any unitary transformation can be built up from combinations 
of the controlled not operation and operations on a single qubit. The implementation 
of a controlled not operation by magnètic resonance involves the preparation of nuclear 
magnetization vectors of a given spin aligned in opposite directions in the transverse plane. 
This "anti-phase" condition, which may be represented in the product operator formalism 
as (say) 2Ij,S2, can be generated in a coupled two-spin system through the initial stages of 



the INEPT pulse sequence |14], relying on (refocused) evolution under the 21^3 z operator 
for a fixed interval l/{2Jjs)- However, the extension of this procedure to more than two 
coupled spins is complicated and not easy to implement. A more direct approach, and 
the one we have employed, is through the use of high-selectivity radio-frequency pulses 
designed to perturb transverse magnetization one line at a time. For example applying a 
vr pulse with Hamiltonian of the form |15| 

+ 2I,R^ + 2S,R^ + 4I,S,R^ (15) 

causes the system to evolve by the unitary operator 

A{2ia^,E,E,E). (16) 

The key observation from the point of view of our work is that more than one such 



line-selective perturbation may be applied simultaneously |16|. Thus any of the unitary 
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operators in Table ^ (and indeed a very wide class of controlled rotations about more gen- 
eral axes) may be produced in the same time that is required to produce the perturbation 



given in (16). It is worth noting that this time is of the same order as that required to 
implement the INEPT sequence. We feel that as weh as being helpful for the present 
work, the method of manipulating spins via simultaneous hne selective pulses may well 
prové advantageous in NMR quantum computers with more spins. 

The experimental task is to shape the radio-frequency pulse envelope so as to achieve suf- 
ficient selectivity in the frequency domain that there is negligible perturbation of the next- 
nearest neighbour of the spin multiplet. In this sense the technique resembles that used 
in pseudo-two-dimensional spectroscopy where the frequency of a soft radio-frequency 
pulse is stepped through the spectrum of interest in very small frequency increments, ex- 
citing the transitions one by one. We investigated several possible pulse shapes for this 
purpose, including rectangular, Gaussian, sine-bell, and triangular, before settling on the 
Gaussian as the most suitable for the task. 

In a weakly-coupled three-spin ISR system the R spectrum is a doublet of doublets with 
splittings Jiü and Jsr- Application of vr pulses to all four transverse i?-spin magneti- 
zation components corresponds to a constant function in the sense of the Deutsch-Jozsa 
algorithm, and the "do nothing" experiment represents the other constant function. The 
balanced functions may be implemented by application of soft vr pulses to the individual 
lines two at a time, for example [0, 0, tt, vr], [0, vr, 0, vr], or [0, vr, vr, 0], where denotes no soft 
pulse. These cases, corresponding to functions /3,/6 and fg have Hamiltonians propor- 
tional to R^. — 2IzTix, R^- — 2S2R;i; and R^ — 4l2S2Ra; respectively. One way to calculate 
the effect of these Hamiltonians is to use Standard product operator manipulations [15U . 
For example one finds that a vr pulse with Hamiltonian of the form R^; — 2Iz^x leaves R^ 
and Sa; unchanged and changes Ix to 21yílx as in Table |. 

The practical implementation is deceptively simple. Starting with a thermal state, a 
hard tt/2 pulse about the y axis (denoted [7r/2]y) excites transverse magnetizations I^;, 
Sx and R^.. The evaluation step is the application of line-selective [njx pulses to the 
individual components of the R multiplet. We may choose to apply soft [vr]^, pulses to all 
four magnetization components, any two of the four, or none at all. In all cases the soft 
pulses are applied simultaneously, while the remaining transitions are simply left to evolve 
freely for the same period of time. However the perturbed magnetization components lose 
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intensity only through spin-spin relaxation during the relatively long interval of the soft 
pulse, T, because the effects of spatial inhomogeneity of the magnètic field are refocused, 
whereas the freely precessing components decay more rapidly, with a shorter time constant 
T2 . This difference in intensities serves to confirm which R transitions were perturbed. 

Experiments were carried out at 400 MHz on a Varian VXR-400 spectrometer equipped 
with a waveform generator which controlled the shaped radiofrequency pulses. The three- 
spin proton system chosen for study was 2,3-dibromopropanoic acid in CDCI3. The three 
coupling constants are Jjr = +11.3 Hz, Jfs = —10.1 Hz, and Jrs = +4.3 Hz. (The 



negative sign of the geminal coupling J/5 |lí] has no particular significance in these exper- 
iments.) Strong coupling effects are evident between spins / and S, with Jis/^is = 0.12. 

Each soft [tï]x pulse can be thought of as acting on one of the four i?-spin magnetizations 
in a rotating frame at the exact resonance frequency of that particular R line. These four 
reference frames rotate at four different freqüències {±Jsr ± Jir)/2 with respect to the 
transmitter frequency centred on the R chemical shift. The x-axes of all four frames must 
be coincident at the beginning of the soft pulse interval T. The duration of the soft pulse, 
T, may be chosen in such a way as to optimize the frequency selectivity. 

The predicted result (Table |5|) , is to convert /- or S'-spin magnetization into various forms 
of multiple-quantum coherence in the six cases where the R magnetization components are 
perturbed in pairs (the balanced functions) but to leave the /- and S'-spin magnetizations 
unaffected in the remaining two cases where the four R magnetization components are all 
perturbed or all left alone (the constant functions). These predictions are clearly borne 
out by the experimental specta shown in the Figure. In principle, complete conversion into 
unobservable multiple-quantum coherence would be detected by the disappearance of the 
appropriate / or S-spin response. In practice, owing to non-idealities of the system (for 
example strong coupling effects between / and S) this is observed as a roughly eightfold 
loss of intensity rather than complete suppression. 

Eight experiments were performed to test the eight cases of Table ^. The transmitter 
frequency was centred on the iï-spin multiplet. Note that the R spectrum remains un- 
perturbed throughout the series, except for the intensity perturbation mentioned above, 
a result of the refocusing effect of the soft vr pulses. The phases of the /- and S signals 
will be determined by the scalar coupling and chemical shift evolution during the period 
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T. These complex phase patterns do not interfere with the Deutsch-Jozsa test because 
this involves only the observation of the "disappearance" of certain signals. These signal 
losses are made clearly evident by displaying absolute-value spectra, which may then be 
integrated. The integrated intensities are shown as percentages of the corresponding in- 
tensities in the top spectrum (no soft-pulse perturbation). Creation of multiple-quantum 
coherence is indicated by the roughly eightfold decrease in intensity in the appropriate 
places; all other I- and S'-spin intensities remain essentially at 100%. This interpretation 
was confirmed in a second experiment with a multiplet-selective soft 7r/2 pulse applied to 
the R spins at the end of the sequence. This has the effect of restoring the "lost" intensities 
by reconverting IR and SR multiple-quantum coherence into observable magnetization. 

Thus, in a single measurement, a distinction can be made between constant and balanced 
functions simply on the grounds of the "disappearance" of /- or S- spin lines. The fact that 
further detalls can be gleaned about the pattern of soft- pulse perturbation is irrelevant to 
the Deutsch-Jozsa algorithm. The extension to systems of more than three coupled spins 
is clear. Because the soft pulses are applied simultaneously, this involves no increase in the 
duration of the perturbation stage. The main limitation would be the magnitude of the 
smallest coupling constant, for this sets the frequency selectivity requirement. Extension 
to more qubits would most likely invoke the introduction of heteronuclear spins such as 
and i^F. 

Acknowledgments The authors are indebted to Dr Èriks Kupce of Varian Associates 
for invaluable advice on the gener ation of shaped selective radio- frequency pulses. 
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Table 1: The eight possible balanced or constant binary functions mapping two bits to 
one bit. 



/ 


u 


/l 


A{E,E,E,E) 


/2 




h 


AiE,E,2a^,2a^) 


h 


A{2a^,2a^,E,E) 


h 


A{2a^,E,2a^,E) 


k 


A{E, 2a^,E,2a^) 


f7 


Ai2a^,E,E,2a^) 


fs 


A{E, 2a^,2aa,,E) 



Table 2: The unitary operators corresponding to the eight constant or balanced binary 
functions mapping two bits to one bit. 
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output 


fl 


+Ix 


+ s. 


— 


+ 2ïx^X 


21a; Ra; 


2Sa;Ra; 


- 4I,S,R, 


f2 


+lx 


+ s. 


— 


+ 2ïxSx 




2Sa;-R'X 


4Ia;Sa;Ra; 


fs 


-Ix 


+ s. 


— R-a; 


2íx^X 




2Sa;R'a; 


+ 4Ia;Sa;Ra; 


/4 


-ïx 




— 


2Ij;Sj; 




2Sa;Ra; 


+ 4Ia:Sa;Ra; 


/5 


+'íx 
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2ïx^x 


- 2I,Rx 


+ 2Sx-RrE 


+ 4Ia:Sa:Ra: 


/6 


+lx 


~ ^x 


— 


2ïx^x 


21a; Ra; 


+ 2Sa:Rx 


+ 4Ia;SxRa; 
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— R-a; 


+ 21a; Sa; 


+ 21a; Ra; 


+ 2Sa;Ra; 


41^ S ^ Jrtí^ 
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— 


+ 21a; Sa; 


~l~ 21a; Ra; 


+ 2Sa;Ra; 


— 4Ia;Sa;Ra; 



Table 3: The output states from a (pseudo) pure initial state after the evaluation stage. 
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Table 4: The effect on input product operators Ix, Sx and R^; of the unitary operators in 
the second column; for example under A{E, 2ax, E, 2ax), the input Sx evolves to 2Sj;Rj; 
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Sx 


Rx 
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Rx 



Table 5: The effect on input product operators Ix, Sx and Rx of the unitary operators in 
the second column; for example under A{E, 2iax,E, 2iax), the input Sx evolves to 2Sy'Rx. 
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Figure 1: Eight absolute-value 400 MHz spectra of 2,3-dibromopropanoic acid obtained 
with the eight different perturbations set out in Table The soft pulses were apphed 
simultaneously with a pulse duration T = 0.65 seconds. Reading from top to bottom, 
these spectra correspond to the functions fi-.-fs of Table |ï[ Integrals of the /- and S-spin 
responses are shown as percentages of those in the top trace. After the evaluation of these 
integrals, the line shapes were improved by pseudo-echo weighting. Note the suppression 
of the appropriate /- and S-spin responses by about an order of magnitude. 
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